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We discuss the possibility of identification of point-like gamma-ray sources (PGS) with small scale dark 
matter (DM) clumps in our Galaxy. Gamma-rays are supposed to originate from annihilation of DM particles 
in the clumps, where annihilation rate is supposed to be enhanced, besides higher density, due to smaller 
relative velocities v of DM particles. We parameterized the annihilation cross section crann('y) in the form of 
an arbitrary power law dependence on the relative velocity v with/without factor of Sommerfeld-Gamow- 
Sakharov, implying existence of a new Coulomb-like interaction. Adopting different parameters of cross 
section and clump, satisfying condition n < 0.2 on density of DM particles of question, they are constrained 
from comparison with Fermi/LAT data on unidentified PGS as well as on diffuse 7-radiation; results are 
applied to concrete DM candidates. Such analysis is found to be sensitive enough to existing uncertainty in 
the density profiles of DM in the clump what can provide a tool for their test. Also we discuss possibilities 
when gamma-radiating clump changes visibly its position on celestial sphere and it is seen as a spatially 
extended gamma-source (EGS), what can be probed in future experiments like Gamma-400. 
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1 Introduction 

From the first articles revealing the indirect effects 
of the cold dark matter (CDM) in the form of 
heavy neutral leptons [1-'^] or supersymmetric par- 
ticles [7, 9, 10], such indirect effects had been the 
subject of intensive studies in the cosmic ray (CR) 
data. It was shown that the DM particles could 
form the hierarchic structures over a wide range of 
scales and masses (from small scale clumps to large 
scale structures) [11-10]. The annihilation of DM 
particles within these structures can give the cosmic 
ray signals [1-8, 17-22]. In the clumps the annihi- 
lation rate should be enhanced due to higher den- 
sity and possibly due to amplification of annihila- 
tion cross section at small relative velocities of DM 
particles, which is especially small in the lightest 
clumps, which are likely to be the most abundant. 
The mentioned factors can lead to that the clumps, 
located in a neighborhood of Solar system, are mani- 
fested as discrete (basically as point-like) gamma-ray 
sources [19,21,23-25]. 



In this paper we continue the previous study of 
given effect [25] and, mainly, make more accurate cal- 
culations, take into account diffuse 7-radiation, con- 
sider possibility of observation of spatially extended 
sources in the light of future experiments. 

2 DM annihilation in the CDM 
clumps 

Predictions of density profile inside the clumps suf- 
fer with some uncertainties [26-29]. For the most 
of calculations in this paper we use profile obtained 
in [19,21,23], which gives rather minimal prediction 
for 7-flux from the clump. Comparison with other 
profile models is given below in terms of the num- 
bers of predicted PGS. The chosen profile (indicated 
further as BGZ) was taken in the form 
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where rj = 1.8, ro = 0.05R, R 10^^{M/Mq)^/^ cm 
with M being the clump mass [19,21]. 

It is supposed that only minor component of 
clumps can survive until the present time [21]. For 
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our estimation ^ = 0.002 is taken as the clumps frac- 
tion from the full density of DM in Galaxy. Here we 
study 7-radiation effect for only minimal clump mass, 
formally assuming that all DM clumps are with this 
mass. In fact, they are predicted to be most abun- 
dant [19,21,30], however effect of clumps (sub-halos) 
from a high-mass tail of mass distribution can be also 
noticeable [ ] . 

A smallness of relative velocity v of DM particles 
concentrated in the clumps may strongly affect to 
annihilation rate if the corresponding cross section 
o'ann depends from v [32,33]. We choose the cross 
section at the parameterized form given below, to 
cover a wide class of models of DM particles: 



O-ann = — 3-C(t;, a). 

v" 



(2) 



Here /3 is a free parameter taking (for generality) the 
continuous range of values (as usual (3 = 1 means 
the s-wave amplitude only, /3 = — 1 means the am- 
plitude from s- and p- waves), (Tq at given /3 is de- 
termined by cosmological density of the particles 0. 
The factor C{v, a) takes explicitly into account a pos- 
sible Coulomb-like interaction (we will refer to it as 
"y-interaction" ) of DM particles, which leads to a 
Sommerfeld-Gamow-Sakharov enhancement [34-37] 
and has the form: 



C {v, a) 
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(3) 



Here a is the fine structure constant of additional in- 
teraction. Such interaction may influence (decrease) 
relic density a little, but may significantly enhance 
the annihilation in the present Universe where the 
particle velocities are small [38-43]. Annihilation ef- 
fects become noticeable even for a subdominant com- 
ponent with Q <C J^cDM as it takes place in case of 
heavy neutrinos [44] to be shown below. 

We suppose that an active (annihilating) compo- 
nent of DM may be both dominant and subdominant, 
i.e. Q < r^cDM with Ocdm ~ 0.2 being the total 
relative density of cold dark matter in Universe. In 
estimation of cosmological density we follow to the 
standard approach [45,46]. It is worth to note that 
the given scheme does not take into account possibil- 
ity of binding pairs of considered particle-antiparticle 
due to y-interaction. The rate of such recombination 
may exceed expansion rate due to high recombination 



cross section, and this process becomes crucial, if DM 
particles decoupled from ambient plasma, leading in- 
evitably to annihilation. It may strongly suppress 
abundance of these DM particles [38]. 

3 DM clumps as gamma-ray 
sources 

In our estimations it is assumed that DM particle of 
question is of Dirac type^ with mass m ~ 100 GeV. If 
annihilating DM particles do not constitute all DM 
{Q < r^cDM) then their contribution to density of 
clumps is assumed to be proportional to il/ficDM 
(eq. (5)). We do not specify annihilation channel 
of photon production, assuming that their averaged 
multiplicity for energy > 100 MeV is A'^^ = 10. 
It is quite typical value for high energy processes at 
respective energy release. 

The photon flux at distance I from the clump center 
is given by 



P 

4^ 



4^ 



(4) 



where the particles/antiparticles number density is 
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Note that the fraction of subdominant DM particles 
should be suppressed in the clumps of mass M < 
Mmin, if they are, where Mmin is the minimal mass 
which could be formed by considered DM particles 
if they prevailed in density. In our study we do not 
take into account this. 

The value ((Tann'w) is determined taking into ac- 
count velocity distribution of DM particles inside the 
clump which is assumed to be Maxwellian one with 
"virial" temperature Tvir = GMm/2R. 

7- Radiation may be registered by LAT, if > 
100 MeV [ ] and their flux exceeds the point source 
sensitivity Fmin ~ 3 • 10~^ cm~^sec~^. The value 
Fjain allows to Calculate the maximal distance /max 
at which the clump can be registered as 7-source. It 
gives for BGZ /max ~ 10~^ pc for /3 = 1.5, ctq = 



^The case of Majorana particles has no differences of prin- 
ciple, except for the case of Coulomb-like interaction which is 
excluded for Majorana particles. 
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Model 


v^UlCU. r 


] Burkert [ ] 


NFW [ ] 


Einasto [ ] Moore [51] 




Isothermal 








1.2 


6.0x10^ 


2.1 X 10^ 


3.3 X 10^ 1.7 X 10^ 



Table 1: The number of PGS as predicted for some density profiles in the unites of that for BGZ profile. 
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^^Mq without y-interaction and 
10 pc with one and the same parameters, 



10 

^max 

and for the Moore model /max ~ 10^^ pc and /max ~ 
10^ pc respectively. All the obtained /max Galactic 
size, what justifies assumption that clump number 
density rid ^ const and corresponds to the local one. 

The number of clumps which may be detected by 
LAT is 

4 n 

Nci = rici • -vr/max> (6) 



where 



^cl 



1.6 • 10' 
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pc 



(7) 



M 

where pioc = 0.3 GeV/cm^. 

The analogous results for some other models are 
given in tab. 1. For all the models "core"-radius and 
clump size were taken the same as in BGZ model. 
It is seen that the BGZ model predicts the minimal 
number of PGSs. So other models become more sen- 
sitive to observational data and the comparison with 
the data becomes very important tool for probing 
them. 

The spatial distribution of clumps at distance /max 
from the Earth is expected to be homogeneous, there- 
fore its distribution on celestial sphere should be 
isotropic. Distribution of unidentified PGS registered 
by Fermi LAT [ ] is almost isotropic except of re- 
gion of the galactic plane. The isotropic component 
include ~100 sources. Supposing that the nature of 
all of them can be related to dark matter clumps, 
it is possible to determine regions of the parameters 
magnitudes /3 and ctq for the typical clump masses 
10"^° IO-'^Mq [22]. The results are shown at fig. 1 
where the factor C is 



m 



^ VlOOGeV/ V 0.2 



CDM 



10 



F„ 



3 • 10 ^ cm ^sec 



This factor includes uncertainties of the chosen pa- 
rameters of DM particle properties (not of clump den- 
sity profile). The predicted number of the clumps 
showing themselves as PGS depends on it as oc 
^-2/3 _ jg ^ggf^i take it into account while con- 
sidering results for other models (see tab. 1). 

The clumps which are not visible (farther than 
/max) should Contribute in diffuse 7-radiation. 7-Flux 
from them in given solid angle can be expressed as 



Fn^rdl 



Pricil 



eff 
halo 



47r 



(9) 



where F and P are introduced in Eq. 4, /halo is the 
distance to the edge of halo along to line of sight and 
^haio ~ effective value (typical for many 

halo density profiles); /max is negligible with respect 
to /halo- One requires that 



^> < $ 



exp 
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(10) 



where <I>exp is the diffuse 7-background measured by 
LAT [5 i]. 

Eq. 10 gives upper limit in the plot fig. 1. As seen, 
the case of clump mass M = 1O~^'^M0 is fully ex- 
cluded, the case of M = 1O~^M0 is constrained but 
up to ~10 PGS are still possible. Higher masses avoid 
such constraining. 

As to the parameters ctq and /?, the range 3T0~^'^ < 
(Jo < 10~^^ cm^ and /3 = 2 (/3 = 1) in case without 
(with) y-interaction are found to be the most inter- 
esting. 

Several specific models of annihilating DM par- 
ticles have been considered: neutralino [55], heavy 
neutrino [44,56] (with y-interaction and without it), 
Kaluza-Klein particles [57], dark atoms OHe [58,59]. 
All candidates except for a heavy neutrino (v^) with 
additional interaction cannot explain LAT data but, 
at the same time, do not contradict them in the case 
of BGZ profile. The case of 1^4 is pointed on fig. 1(b). 
As seen, the part of nonidentified point-like LAT 



3 




(a) The case without y-interaction is shown. (b) The case with y-interaction is shown. 



Figure 1: Allowed and forbidden regions of parameters /3 and ctq are shown as they obtained on the base of 
Fermi LAT data on PGS and diffuse 7-radiation. Two values of the minimal clump mass are given. 1^4 on 
fig. 1(b) corresponds to the cases of heavy neutrino in the mass range 45-49 GeV. 



sources may be explained due to z/4p4-annihilation at 
mass niui ~ 46 — 49 GeV. Typical 7-spectrum from 
47 GeV neutrinos annihilation is shown at the fig. 2 
in comparison with measured spectrum of one of the 
non-identified PGS (annihilation spectrum was ob- 
tained with the help of Monte-Carlo generator Pythia 
6.4 [60]). It is worth to note that the heavy neutrino 
parameters are strongly restricted by underground 
experiments [44,61]. At the same time the predic- 
tion of z^4 relic density suffers with uncertainty con- 
nected with recombination of the pair of y-interacting 
neutrino-antineutrino in early Universe [44]. 

4 Some possible features of the 
model 

4.1 Spatially extended sources 

The distribution of clumps does not exclude a possi- 
bility of their presence in close vicinity to Solar sys- 
tem. Therefore an angular size of such the sources 
can exceed the angular resolution of detector 5 (cor- 
responding to the solid angle, say, = 7r((5/2)^). 

In the capacity of criterion of extended gamma- 
source (EGS) we choose the following. The minimal 
amount of photons needed to recognize a source over 
background from a region is defined by the min- 
imal point source flux -Fmin- Statistically significant 
signal from region with the size of I'Qs must be then 




E, GeV 

Figure 2: Expected spectrum from i/41/4-annihilation 
is shown in comparison with the observed spectrum 
of nonidentified source 2FGL J1653. 6-0159. 

at least > y/i'Fram- For clump candidate to EGS we 
require fulfilment of similar condition for the circle 
region of vrJ^ = 4^2^ size around the clump center 
excluding the 7r(5/2)^ center part: 

Fi^ns) - F{ns) > 2Fn,in. (11) 

Here for minimal flux a little harder condition was 
taken than it should be proceeding from the num- 
ber of ils regions which the ring [6/2 .. .5] covers, i.e. 
\/3Fmin. However, this criterion (11) is softer than 
requirement ^(^^,5) > i^min for several adjacent re- 
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M/Mq 

Figure 3: Allowed and forbidden regions of param- 
eters P and M are shown for two typical values of 
ctq without y-interaction. The solid and dashed near 
red lines show 1 extended source for Fermi LAT and 
Gamma-400 detectors respectively. 

gions simultaneously or even F{vVts) ^ T^Fmin for all 
the uVts region wholly, and it needs a special analysis 
of observation data. 

Based on Eq. (11) we find that Fermi LAT can 
hardly observe any clump as non-point-like source, 
what does agree with observation data. However, it 
becomes possible at the improved in future experi- 
ment values of -Fmin and 5. At the fig. 3 we show 
the regions of parameters of clump mass M and (3 
satisfying to data on diffuse and discrete sources of 
7-radiation. Dashed line shows EGS discovery po- 
tential at the achievement oi 5 = (5g4oo ~ 0.3° (vs 
(5lat ~ 0.6° for Fermi LAT at > 1 GeV [52]) and 
the same -Fmin as for Fermi LAT. In principle, such 
values can be reached for Gamma-400 detector [62] 
at some energy interval. EGS discovery will help to 
test the clump model. 

Note, that improvement (minimization) of i^min (in 
perfect case, minimal flux per unit of solid angle) is 
found to be more promising from viewpoint of EGS 
discovery than improvement of 6 only. However, at 
better 6 one can partially compensate low sensitivity 
to photon flux by taking several fJ^ region. 

4.2 "Traveling" sources 

Provided DM clump accounts for PGS, there can 
be sources which change visibly its position on ce- 



lestial sphere within time of observations (the time 
between EGRET experiment [()•"')] and Fermi LAT 
is ~10 years), due to its closeness to Solar sys- 
tem. Clump velocity in Galaxy is f ~ 250 km/s. 
For 10 years they travel the distance about A ~ 
2.6 • 10^^ pc. Let us estimate the number of clumps 
with visible movement. The number of clumps 
per the spherical layer r ^ r + dr and per interval 
cos^ cos^ -|- dcosO, where 9 is the angle between 
radius- vector r (originating from observer) and vec- 
tor of clump velocity v, is 



dN 



drd cos ( 



47rr^nci-. 



(12) 



Here we refer to angular distribution of v as isotropic 
one with ranging from to tt. 

Integration (12) gives the distribution in angle 



Asin6' 



which clump travels on the sky. 



dN 



dN 



drd cos 6 



drd cos 95 (</) — 0(r, 6)) 



3 



(13) 



Here we changed the variables r and cos 6 to (j) 
through integration with 5-function (divergence at 
i?i> — 7- is a consequence of our approximation rif-i = 
const for < r < oo). 

Estimation of the number of the moved (visibly) 
sources gives 



N, 



dN 



1 



A3 



< 0.4 



< 



IO-'^Mq 



M 



(14) 



'^EGRET + ^LAT ~ 



Here 0max ^ 1 does not afi^ect the result, 
max((5i.esoi,5iim), where 

^resol ~ 

1.47° with Regret » 1-34° and (5lat » 0.6° be- 
ing the angular resolutions of detectors EGRET and 
LAT for ^ 1 GeV, ~ A/^max corresponds 
to angular shift of the clump at maximal observ- 
able distance. For maximal estimations in (14) we 
put 4>rmn = ^resoh that gave result independent on 
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the model of the clump density profile. In case of 
<?^'mm = <5iim > '^rcsoh the visible clumps (if they 
are) should be traveling. 

Thus, as seen from (14), for mass of clumps M ~ 
1O~^M0 the probability to find the moved sources be- 
comes noticeable (for any density profile). The value 
M ^ Mq is close to the minimal one which is 
allowed from data on diffuse 7-background (fig. 3). 
However, time of observation(s) grows and future ex- 
periment like Gamma-400 will work at better angular 
resolution. Both from the time and from the resolu- 
tion the expected number of traveling sources (TS) 
depends as a cube power (see (14)). It challenges for 
future experiment to probe the model with searching 
for TS for clump mass up to ~ (10~^ H- 1Q~^)Mq. 
Even now, one may consider time interval of 20 years 
(from the start of EGRET) and expect to find several 
TS. 

In case of A'trav ^ 1 the movement of the sources 
should have a regular character accounted for by 
solar system motion around Galactic center (GC). 
Respective data analysis would provide much more 
strong test for discussed model, since virtually all the 
known astrophysical explanations of PGSs do not re- 
late them with objects of halo (but rather with galac- 
tic disc or distant galaxies, in case of which observa- 
tions are insensitive to solar system motion around 
GC). 

Existing difference in data on non-identified PGSs 
of EGRET and LAT (LAT confirms existence of only 
30-40% of the sources from catalogue 3EG [53, 63] 
and observes new sources which had not been reg- 
istered by EGRET) could be explained by effects of 
traveling clump-sources in small part. 

5 Conclusion 

In this paper we refine the previous results [ -":] and 
have shown that DM clumps could be point-like (and 
extended) sources of the 7-radiation and they can 
partially explain non-identified 7-sources, registered 
by LAT and EGRET. The proposed method allows 
to estimate an amount of accessible to observation of 
DM clumps for various models of DM particle and 
density profile. Note that the suppression of the sub- 
dominant fraction of DM particles in clumps of mass 
M < Mmin (if they are) has not been taken into ac- 
count and requires special research. 



The values of parameters ctq and /?, at which DM 
model is either consistent or inconsistent with the 
LAT data, have been determined for the most "con- 
servative" model BGZ [21]. The high sensitivity of 
the predictions to the choice of density profile model 
is shown. 

The clumps, situated in a close vicinity of Solar 
system, may account partially for noncoincidences 
between catalogues EGRET and LAT (the sources, 
registered by EGRET and not confirmed by LAT) 
and also for a spatially extended gamma-ray sources 
which can be detected by future gamma-telescopes. 

Possibilities that gamma-radiating clump changes 
visibly its position on celestial sphere and it is seen 
as a spatially extended gamma-source can play a role 
of distinctive features, to be probed in future experi- 
ments. It would allow to distinguish from alternative 
models of possible gamma-source origin and in par- 
ticulary the model associated with primordial black 
hole clusters [64,65]. 
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